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ABSTRACT 

We investigate morphological structure parameters and local environments of distant moderate- 
luminosity active galactic nucleus (AGN) host galaxies in the overlap between the HST/ ACS ob- 
servations of the Great Observatories Origins Deep Survey (GOODS) and the two Chandra Deep 
Fields. We compute near-neighbor counts and BViz asymmetry [A) and concentration (C) indices 
for K, 35,500 GOODS/ACS galaxies complete to zgso ~ 26.6, including the resolved hosts of 322 
X-ray-selected AGNs. Distributions of (1) zsso asymmetry for 130 zgso < 23 AGN hosts and (2) 
near- neighbor counts for 173 zgso < 24 AGN hosts are both consistent with non-AGN control sam- 
ples. This implies no close connection between recent galaxy mergers and moderate-luminosity AGN 
activity out to appreciable look-back times (z < 1.3), approaching the epoch of peak AGN activity 
in the universe. The distribution of zg5o C for the AGN hosts is offset by AC w -fO.5 compared 
to the non-AGN, a 6.4 a discrepancy much larger than can be explained by the possible influence 
of unresolved emission from the AGN or a circumnuclear starburst. The local universe association 
between AGN and bulge-dominated galaxies thus persists to substantial look-back time. We discuss 
implications in the context of the low-redshift supermassive central black hole mass correlation with 
host galaxy properties, including concentration. 

Subject headings: galaxies: active — galaxies: structure — surveys — X-rays: galaxies 



1. INTRODUCTION 

The connection between active galactic nuclei (AGNs) 
and their host galaxies, and the evolution in that re- 
lationship over cosmic time, have attracted great inter- 
est in recent years. This includes the discoveries that 
most nearby massive galaxies harbor central supermas- 
sive black holes (SMBHs; Magorrian et al. 1998), that 
AGNs in the local universe (z < 1) reside predominantly 
in massive, bulge-dominated host galaxies (Kauffmann 
et al. 2003), and that a tight correlation exists locally 
between SMBH mass and host galaxy properties such as 
bulge velocity dispersion and light-profile concentration 
(Ferrarese & Merritt 2000; Gebhardt et al. 2000; Gra- 
ham et al. 2001). The 1-2 Ms Chandra Deep Fields 
(CDF-N and CDF-S; Brandt et al. 2001; Giacconi et 
al. 2002), which have now resolved much of the cosmic X- 
ray background into moderate- luminosity AGNs a,t z \ 
(Alexander et al. 2003; Barger et al. 2003; Szokoly et 
al. 2004), provide a unique AGN sample to probe these 
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locally observed SMBH-host relationships out to epochs 
nearing the peaks of star formation and AGN activity 
in the universe. This is one of the aims of the Great 
Observatories Origins Deep Survey (GOODS; Giavalisco 
et al. 2004), which has obtained deep multicolor Hubble 
Space Telescope {HST) Advanced Camera for Surveys 
(ACS) image mosaics across the most sensitive regions of 
the CDF areas. In this Letter we report on the local en- 
vironments and rest-frame optical morphologies of AGN 
host galaxies in the GOODS fields as well as the implica- 
tions for SMBH-galaxy coevolution and the merger- AGN 
connection. 

The largest pre-GOODS investigation of HST-imaged 
CDF sources was Koekemoer et al. (2002), with 41 CDF- 
S 1 Ms sources in three moderately deep Wide Field 
Planetary Camera 2 (WFPC2) pointings. Grogin et 
al. (2003, hereafter G03) studied the HST morphologies 
and local environments of these faintest X-ray sources as 
compared with the X-ray undetected population. The 
AGNs were preferentially located in galaxies with highly 
concentrated light profiles, but the AGN hosts could 
not be differentiated from the non-AGN based on light- 
profile asymmetry or frequency of near-neighbors. The 
G03 conclusions were as follows: (1) distant moderate- 
luminosity AGNs did not show a connection between 
recent (< 1 Gyr) galaxy merger/interaction and AGN 
activity; and (2) the z ~ 1 galaxy population already 
showed evidence of the SMBH-bulge correlation. Now 
that deeper and much larger area GOODS imaging is 
available across both CDFs, we verify these results with 
much larger samples of both CDF AGNs and quiescent 
galaxies. We also place new constraints on the evolution 
in merger- AGN connection and SMBH-bulge correlation 
with the extensive redshift information now accumulated 
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in these fields. We adopt a cosmology with Hq = 70 km 
Mpc-\ fl,n = 0.3, and = 1 - = 0.7. Magni- 
tudes are given in the AB system. 

2. OBSERVATIONS AND SAMPLE SELECTION 

Our analyses employ the HST/ ACS image mosaics in 
F606W {V), F775W (i), and F850LP (zgso) from the first 
three epochs of GOODS HST observations in both the 
northern ("GOODS-N") and southern ("GOODS-S") 
fields (Giavalisco et al. 2004). We use the zgso-detected 
source catalog detailed in Ravindranath et al. (2004), 
trimmed by applying minimum thresholds in compact- 
ness (> 4 pixels) and signal-to-noise ratio (> 5) opti- 
mized for removing spurious detections. The remaining 
16,632 GOODS-S sources and 18,878 GOODS-N sources 
are complete to zgso ~ 26.6 and form the basis of our 
environmental and structural parameter analyses. 

We identify candidate AGN hosts from source cata- 
logs of the 1 Ms CDF-S and 2 Ms CDF-N X-ray im- 
ages reduced and source-extracted in a consistent fashion 
(Alexander et al. 2003). These two X-ray surveys provide 
the deepest views of the universe in the 0.5-8.0 keV band. 
Within the respective GOODS-N(S) areas, the CDF- 
N(S) contains 324(223) X-ray sources down to compara- 
ble sensitivity limits of w 1.0(1.3) x 10^^® ergs s 
at 0.5-2.0 keV and « 7.2(8.9) x lO-^^ 
at 2-8 keV. Coordinate-matching to the zgso catalog 
yields unambiguous counterparts for > 80% of the CDF 
sources (F. E. Bauer et al. 2005, in preparation). Many 
are comparatively nearby and optically bright starbursts 
and "quiescent" galaxies contaminating our desired X- 
ray-selected AGN sample. The extensive redshift cov- 
erage of ^850 ^ 24 CDF counterparts allows us to ex- 
clude these non-AGNs with a luminosity threshold of 
L(2-8 keV) > 10'*^ ergs s~^ The resuhing Lx-limited 
AGN sample of 322 galaxies contains few CDF sources at 
z < 0.4, so we choose this as our lower limit for redshift- 
evolution analyses (f|5)|. 

We investigate evolutionary trends in morphology and 
environment among the GOODS AGN and non-AGN 
populations by constructing complete volume-limited 
subsamples. We estimate absolute magnitudes for the 
GOODS-S galaxies by using the photometric redshift 
database of Mobasher et al. (2004), who claim an ac- 
curacy of Az/(1 + z) < 0.1 at Z850 < 24.5 for AGN and 
non-AGN alike. Although we lack a comparable pho- 
tometric redshift database for the GOODS-N field, red- 
shift measurements of the CDF-N sources (Fernandez- 
Soto, Lanzetta, & Yahil 1999; Cohen et al. 2000; Barger 
et al. 2003) are largely complete to a similar depth 
(•Z850 ^ 24.5). Hence, we include another 135 CDF-N 
sources with measured redshifts for our volume-limited 
AGN sample. 

To probe the evolution of the optical morphology of 
AGN hosts versus the field, we compute rest-frame B- 
band (hereafter Bq) quantities out to the limit z < 1.3 
accessible to our reddest filter (zsso)- To estimate Mb 
and the Bq structural parameters, we linearly interpo- 
late between V and i for galaxies at 0.31 < z < 0.73, 
between i and zgso for 0.73 < z < 1.09, and use zgso 
quantities for 1.09 < z < 1.3. The redshift survey limit 
of Z850 « 24.5 corresponds to Mb ~ —19.5 at z = 1.3. 
We adopt this as our limiting absolute magnitude, satis- 
fied by 1090 GOODS-S galaxies in the non-AGN sample 



and another 37(42) AGN hosts from CDF-N(S). These 
volume-limited samples probe to w 1 mag fainter than 
and thus arc not restricted to the highest luminosity 
galaxies. Coincidently, z 1.3 is the limit for Chandra 
detection of L{2-8 keV) > 10'*^ ergs s^^ sources through- 
out the GOODS regions. As a result, our AGN sample 
is essentially complete. 



3. CONCENTRATION AND ASYMMETRY INDICES 

We quantify the GOODS galaxy morphologies via non- 
parametric indices of concentration C and asymmetry 
A (Conselice 2003 and references therein). Index C 
scales with the ratio of radii containing 80% and 20% 
of a source's total flux, C = 5 log(ro.8/?'o.2), and in- 
creases toward bulge-dominated morphologies. Index 
A is the flux-normalized residual of the source pixels 
S differenced with their 180°-rotated counterpart 5*180 : 

A = min (j2pi^ \S - Siso\ Epix l^l) ~ ^o, where Aq is 
a background term and the minimization is over a 0.2 
pixel grid of possible centers of rotation. While A mod- 
erately increases toward disk-dominated morphologies, it 
is driven to large values {A > 0.35) for galaxies with re- 
cent or ongoing interaction. 

Figure ^ shows the zgso indices C (top panel) 
and A {bottom panel) for resolved sources in 
both GOODS fields versus zsso magnitude. 



O 

' I — I 



o 
o 



0.8 

£"0.6 
■+-> 

CD 

gO.4 
^ 



-0.2 




□CDF-N 



H — I h 



-I — h 



1 



■ACDF-S 




18 



SO 



22 



24 



26 



z,,, [ABMAG] 



Fig. 1. — Concentration index (top) and asymmetry index 
(bottom) vs. magnitude as measured in zg^g for resolved GOODS 
sources. The large symbols represent the AGN sample in the north 
{red squares) and south (blue triangles); the small dots represent 
the non-AGNs in the north (yellow) and south (magenta). The 
connected green crosses and their error bars denote the median 
values and measurement errors for the non-AGNs in successive 1 
mag bins. 
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The asymmetries of the AGN hosts {large symbols) and 
non-AGN {small dots) are not clearly separated. How- 
ever the AGN host concentrations clearly populate the 
upper end of the field distribution throughout the regime 
of good signal-to- noise ratio (zgso < 24). Such large C- 
values are preferentially associated with massive early- 
type galaxies, although we note that the AGN hosts 
span a broad range of morphology (e.g., Koekemoer et 
al. 2002). The declining C with host magnitude is in- 
terpreted as an evolutionary effect; galaxies at higher 
redshift (and thus generally fainter) are intrinsically less 
concentrated (Conselice et al. 2003). When compar- 
ing the C distributions for zgso-limited samples of AGN 
and non-AGN f ti5.1|l . we compensate for this magnitude- 
dependent bias in C by resampling the non-AGN popu- 
lation to match the AGN hosts' magnitude distribution. 



4. NEAR-NEIGHBOR FREQUENCY 

In assessing the role of environment in AGN activity, 
we complement the A analysis with a comparison of the 
near-neighbor counts around zgso < 24 AGN hosts versus 
the non-AGN. Qualifying neighbors must satisfy both 
proximity and relative brightness criteria. We investigate 
two different definitions of proximity threshold d: one 
that scales with the Petrosian radius rp of the primary 
galaxy, d < 3rp, and another that remains fixed for all 
galaxies, d < 8". The latter choice, consistent with the 
analysis of G03, corresponds to 54 kpc a,t z = 0.6 and 
varies by only zt25% over the range 0.4 < z < 1.3. 

To limit contamination of the neighbor statistics by 
chance superpositions of background galaxies, a neigh- 
bor is rejected if more than 2 mag fainter than the com- 
parison galaxy. This relative magnitude cutoff, in the 
presence of steeply increasing galaxy number counts with 
magnitude, introduces a bias towards more neighbors 
around fainter galaxies. When comparing AGN and non- 
AGN near-neighbor counts, we therefore resample the 
non-AGN to match the magnitude distribution of the 
AGN hosts, analogous to the procedure used in compar- 
ing C distributions (fj^j). Because the resolved CDF opti- 
cal counterparts typically show only minor ffux contribu- 
tion from the active nucleus (G03), the relative faintness 
of qualifying AGN neighbors is not significantly biased 
with respect to the non-AGN. We discuss the similar- 
ity of the AGN and non-AGN near-neighbor frequency 
histograms in ti5.2l 



5. DISCUSSION 

5.1. Bulge-dominated AGN Hosts at z ^ 0.4-1.3 

The top left panel of Figure |21 notes the Kolmogorov- 
Smirnov (K-S) test probabilities for the null hypothesis 
that the zgso < 23 AGN host and non-AGN C-values 
could be drawn from the same underlying distribution. 
The quoted probabilities are the median values from 
1000 resamplings of the non-AGN sample matched to 
the magnitude distribution of z < 23 AGN hosts to 
remove magnitude-dependent C bias (see iJ^J- The C 
distributions are highly inconsistent at the 6.4 a level 
(Pk-s = 1-6 X 10^^"). Moreover, both the northern 
and the southern zgso < 23 AGN host concentrations are 
individually discrepant with their respective non-AGN 
counterparts at Pk-s < 10^^. 
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Fig. 2. — Cumulative distribution functions of concentration 
index (left panels) and asymmetry index (right panels) for resolved 
GOODS AGN hosts from the CDF-S (dotted blue line), CDF-N 
(dashed red line), and the combination of both CDFs (thick solid 
black line) for a flux-limited sample (top row) and a volume-limited 
sample divided into three redshift bins (bottom three rows). The 
indices are compared in zgso for the flux-limited sample, and in 
(interpolated) rest-frame B for the volume-limited sample. Null- 
hypothesis probabilities P from K-S comparisons with the control 
sample of non-AGN GOODS-South galaxies (thin solid black line) 
are also noted in each panel. 

The median C offset of -1-0.5 is consistent with the 
G03 measurement based on / < 23 sources in three 
ffST/WFPC2 pointings in the CDF-S (including 25 X- 
ray-detected sources) and far exceeds the '^0.1 enhance- 
ment in C expected from nuclear point-source optical 
fiux in the AGN sample (see G03). Moreover, C en- 
hancement by potential circumnuclear starbursts is dis- 
counted because (1) unresolved starbursts reduce to the 
previous case for AGN point-source contamination and 
(2) resolved (> 500 pc) central starbursts with sufficient 
fiux to bias C would be inconsistent with the low in- 
cidence of starburst-type spectral energy distributions 
(SEDs) among the GOODS AGN hosts as compared to 
the field (Mobasher et al. 2004). 

The traditional conception of "local AGN = Seyfert 
= late-type host" has been refuted by the Kauffmann 
et al. (2003) analysis of thousands of spectroscopically 
identified low-rcdshift {z < 0.3) AGN host galaxies from 
the Sloan Digital Sky Survey (SDSS). Local AGN of all 
luminosities reside almost exclusively in massive hosts 
with sizes, stellar mass densities, and concentration in- 
dices similar to ordinary early-type SDSS galaxies. The 
enhanced C among the GOODS AGN hosts now indi- 
cates that nuclear activity remains preferentially asso- 
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dated with bulge-dominated galaxies out to substan- 
tial lookback times {z < 1.3). Our z ^ 0.4-1.3 sam- 
ple largely bridges the span of cosmic time between the 
quasar epoch, where accretion-driven luminosity is dom- 
inated by high-mass SMBHs, and the recent-epoch AGN 
hosts probed by SDSS. If the locally-observed tight cor- 
relation between SMBH mass and host C (Graham et 
al. 2001) similarly extends to z ~ 0.4-1.3, then our re- 
sults newly suggest that the accretion-driven luminosity 
of the universe is dominated by the most massive SMBHs 
at virtually all times. 

The discovery of an epoch beyond which SMBH mass 
and host-galaxy properties (including C) lose their tight 
correlation could strongly constrain theories of SMBH 
and host-galaxy co-evolution. To ascertain if the C dis- 
crepancy between AGN and non-AGN shows any evo- 
lution with redshift, we divide our volume-limited sam- 
ple into three redshift bins spanning « 550 Mpc'^ each: 
0.4 < z < 0.85, 0.85 < z < 1.1, and 1.1 < z < 1.3. 
The AGN host C is clearly elevated in all three bins 
(Fig. 121 three bottom left panels), reflected in the persis- 
tently low logPK-s 6.8 to -4.7. Thus our GOODS 

AGN host-galaxy sample populates the high end of both 
the Lx distribution (by construction) and the C distri- 
bution throughout the range 0.4 < z < 1.3. Pushing this 
analysis beyond z « 1.3 faces multiple obstacles, includ- 
ing (1) identification of obscuration-unbiased moderate- 
luminosity AGNs {Lx > 10^^ ergs s^^), requiring Chan- 
dra exposure depths of several megaseconds; (2) con- 
struction of a large, complete Mb < —19.5 field sample 
in the so-called redshift desert; and (3) well-resolved rest- 
frame optical light profiles out to meaningful isophotes, 
requiring ^ O'.'l resolution J-band imaging at extreme 
depths to overcome the sharply increasing surface bright- 
ness dimming. 

Although tantalizing to conclude that the Graham 
et al. (2001) SMBH-bulge correlation now persists to 
z < 1.3, it is unclear whether Lx is a reasonable proxy 
for SMBH mass at these rcdshifts. At low redshift, AGNs 
with well-constrained SMBH mass are observed to have 
(1) Lx ^ ^Edd oc Mbh, but with a large scatter to lower 
Eddington ratios (Woo & Urry 2002) and (2) an AGN 
fundamental plane in Lx, L Radio, and Mbh, but poor 
correlation between Lx and Mbh individually (Merloni 
et al. 2003). However Barger et al. (2005) have recently 
claimed a tight L^-Mbh correlation at z ~ 1. Direct 
measurement of the SMBH masses of GOODS galaxies 
would be ideal but extremely challenging for moderate- 
luminosity (and often dust-obscured) AGNs at these red- 
shifts. 

5.2. No Merger- AGN Connection at z ^ 0.4-1.3? 

Unlike the disparate C distributions, the asymmetry 
indices of the 130 resolved zgso < 23 GOODS AGN 
hosts are statistically indistinguishable from the non- 
AGN (Pk-s = 0.97). This result reinforces the G03 
findings based on a subset of 25 X-ray-detected sources 
with / < 23. Figure |21 {right panels) shows that this 
similarity in A exists in Bq throughout all three redshift 
bins and separately among northern and southern AGN 
subsamples. The lowest of the K-S test probabilities, for 
the combined AGN sample at 0.4 < z < 0.85, does not 
exceed a 2 cr rejection of the null hypothesis. 

Recent/ongoing galaxy mergers in the local universe 



generally have large A enhancement {A > 0.35; Gon- 
selice et al. 2003). Furthermore, A^-body simulations 
suggest that even minor mergers have significant A en- 
hancement up to 1 Gyr from the onset of the merger 
(Walker, Mihos & Hernquist 1996). The lack of differ- 
entiation between AGN and non-AGN A distributions 
therefore implies that recent merging/interaction is no 
more prevalent among the AGNs. This in turn argues 
against the hypothesis that AGN fueling is directly linked 
to recent (< 1 Gyr) galaxy merging/interaction, while fa- 
voring mechanisms such as low-level gas accretion from 
the intergalactic medium and/or bar instability in disks. 

HST imaging of 20 nearby (z < 0.3) quasars by Bahcall 
et al. (1997) suggested that galaxy mergers/interactions 
are relatively common among the highest luminosity 
AGNs (>10''^'^ ergs s~^). However, our finding of no A 
enhancement among more distant (z ^ 0.4-1.3) AGNs 
at lower luminosities (^lO'*'^'^ ergs s^^) has precedent in 
the comparable low-redshift sample analyzed by Corbin 
(2000). The GOODS combination of deep Chandra and 
HST imaging now suggests a persistent merger- AGN dis- 
connect among moderate-luminosity AGNs out to look- 
back times nearing the peak of AGN activity in the uni- 
verse. The limited GOODS solid angle provides too few 
high-luminosity AGNs to test the Bahcall et al. (1997) 
conclusions over the same redshift range. This may be 
remedied by the wider area Galaxy Evolution from Mor- 
phology and SEDs project HST imaging within the ex- 
tended GDF-S, where Sanchez et al. (2004) have already 
noted a higher incidence of merger/interaction among 15 
optically selected, z ~ 0.5-1.1, type 1 AGNs with lumi- 
nosities spanning the Seyfert/quasar boundary. 

Our comparison of the near-neighbor frequency his- 
tograms for AGNs and non-AGNs yields a large x^-test 
probability of the null hypothesis both for the d < 3rp 
threshold [P(x^) = 0.84] and the d < 8" threshold 
[P{x'^) — 0.58]. We note that the ~ 2(t discrepancy in 
near-neighbor frequency previously observed in G03 is no 
longer apparent for the substantially enlarged samples of 
the current work. We conclude that local environment, 
like host asymmetry, is not well correlated with AGN ac- 
tivity. This result, among moderate-redshift AGNs, now 
extends similar findings of environment- AGN disconnect 
at low redshift from recent analyses of the Southern Sky 
Redshift Survey (Maia, Machado, & Willmer 2003) and 
the SDSS (Miller et al. 2003). 
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